Search for: All records

Antibiotic resistance (AR) presents a global health challenge, necessitating an improved understanding of the ecology, evolution, and dissemination of antibiotic resistance genes (ARGs). Several tools, databases, and algorithms are now available to facilitate the identification of ARGs in metagenomic sequencing data; however, direct annotation of short-read data provides limited contextual information. Knowledge of whether an ARG is carried in the chromosome or on a specific mobile genetic element (MGE) is critical to understanding mobility, persistence, and potential for co-selection. Here we developed ARGContextProfiler, a pipeline designed to extract and visualize ARG genomic contexts. By leveraging the assembly graph for genomic neighborhood extraction and validating contexts through read mapping, ARGContextProfiler minimizes chimeric errors that are a common artifact of assembly outputs. Testing on real, synthetic, and semi-synthetic data, including long-read sequencing data from environmental samples, demonstrated that ARGContextProfiler offers superior accuracy, precision, and sensitivity compared to conventional assembly-based methods. ARGContextProfiler thus provides a powerful tool for uncovering the genomic context of ARGs in metagenomic sequencing data, which can be of value to both fundamental and applied research aimed at understanding and stemming the spread of AR. The source code of ARGContextProfiler is publicly available atGitHub. 

Despite the recent surge of viral metagenomic studies, it remains a significant challenge to recover complete virus genomes from metagenomic data. The majority of viral contigs generated from de novo assembly programs are highly fragmented, presenting significant challenges to downstream analysis and inference. To address this issue, we have developed Virseqimprover, a computational pipeline that can extend assembled contigs to complete or nearly complete genomes while maintaining extension quality. Virseqimprover first examines whether there is any chimeric sequence based on read coverage, breaks the sequence into segments if there is, then extends the longest segment with uniform depth of coverage, and repeats these procedures until the sequence cannot be extended. Finally, Virseqimprover annotates the gene content of the resulting sequence. Results show that Virseqimprover has good performances on correcting and extending viral contigs to their full lengths, hence can be a useful tool to improve the completeness and minimize the assembly errors of viral contigs. Both a web server and a conda package for Virseqimprover are provided to the research community free of charge. 

ML Framework for PPCPs fate in WWTPs. 

We study in detail the influence of different chemical potentials (baryon, electric charge, strange, and neutrino) on how and how fast a free gas of quarks in the zero-temperature limit reaches the conformal limit. We discuss the influence of non-zero masses, the inclusion of leptons, and different constraints, such as charge neutrality, zero-net strangeness, and fixed lepton fraction. We also investigate for the first time how the symmetry energy of the system under some of these conditions approaches the conformal limit. We find that the inclusion of all quark masses (even the light ones) can produce different results depending on the chemical potential values or constraints assumed. A positive or negative deviation of 10% from the pressure of free massless quarks with the same chemical potential was found to take place as low as μB=77 to as high as 48,897 MeV. This illustrates the fact that the “free” or conformal limit is not a unique description. Finally, we briefly discuss what kind of corrections are expected from perturbative QCD as one goes away from the conformal limit. 

Cu is an antimicrobial that is commonly applied to premise (i.e., building) plumbing systems for Legionella control, but the precise mechanisms of inactivation are not well defined. Here, we applied a suite of viability assays and mass spectrometry-based proteomics to assess the mechanistic effects of Cu on L. pneumophila. Although a five- to six-log reduction in culturability was observed with 5 mg/L Cu2+ exposure, cell membrane integrity only indicated a <50% reduction. Whole-cell proteomic analysis revealed that AhpD, a protein related to oxidative stress, was elevated in Cu-exposed Legionella relative to culturable cells. Other proteins related to cell membrane synthesis and motility were also higher for the Cu-exposed cells relative to controls without Cu. While the proteins related to primary metabolism decreased for the Cu-exposed cells, no significant differences in the abundance of proteins related to virulence or infectivity were found, which was consistent with the ability of VBNC cells to cause infections. Whereas the cell-membrane integrity assay provided an upper-bound measurement of viability, an amoebae co-culture assay provided a lower-bound limit. The findings have important implications for assessing Legionella risk following its exposure to copper in engineered water systems. 

Abstract With growing calls for increased surveillance of antibiotic resistance as an escalating global health threat, improved bioinformatic tools are needed for tracking antibiotic resistance genes (ARGs) across One Health domains. Most studies to date profile ARGs using sequence homology, but such approaches provide limited information about the broader context or function of the ARG in bacterial genomes. Here we introduce a new pipeline for identifying ARGs in genomic data that employs machine learning analysis of Protein-Protein Interaction Networks (PPINs) as a means to improve predictions of ARGs while also providing vital information about the context, such as gene mobility. A random forest model was trained to effectively differentiate between ARGs and nonARGs and was validated using the PPINs of ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, andEnterobacter cloacae), which represent urgent threats to human health because they tend to be multi-antibiotic resistant. The pipeline exhibited robustness in discriminating ARGs from nonARGs, achieving an average area under the precision-recall curve of 88%. We further identified that the neighbors of ARGs, i.e., genes connected to ARGs by only one edge, were disproportionately associated with mobile genetic elements, which is consistent with the understanding that ARGs tend to be mobile compared to randomly sampled genes in the PPINs. This pipeline showcases the utility of PPINs in discerning distinctive characteristics of ARGs within a broader genomic context and in differentiating ARGs from nonARGs through network-based attributes and interaction patterns. The code for running the pipeline is publicly available athttps://github.com/NazifaMoumi/PPI-ARG-ESKAPE 

Abstract BackgroundWhile there is increasing recognition of numerous environmental contributions to the spread of antibiotic resistance, quantifying the relative contributions of various sources remains a fundamental challenge. Similarly, there is a need to differentiate acute human health risks corresponding to exposure to a given environment, versus broader ecological risk of evolution and spread of antibiotic resistance genes (ARGs) across microbial taxa. Recent studies have proposed various methods of harnessing the rich information housed by metagenomic data for achieving such aims. Here, we introduce MetaCompare 2.0, which improves upon the original MetaCompare pipeline by differentiating indicators of human health resistome risk (i.e., potential for human pathogens to acquire ARGs) from ecological resistome risk (i.e., overall mobility of ARGs across a given microbiome). ResultsTo demonstrate the sensitivity of the MetaCompare 2.0 pipeline, we analyzed publicly available metagenomes representing a broad array of environments, including wastewater, surface water, soil, sediment, and human gut. We also assessed the effect of sequence assembly methods on the risk scores. We further evaluated the robustness of the pipeline to sequencing depth, contig count, and metagenomic library coverage bias through comparative analysis of a range of subsamples extracted from a set of deeply sequenced wastewater metagenomes. The analysis utilizing samples from different environments demonstrated that MetaCompare 2.0 consistently produces lower risk scores for environments with little human influence and higher risk scores for human contaminated environments affected by pollution or other stressors. We found that the ranks of risk scores were not measurably affected by different assemblers employed. The Meta-Compare 2.0 risk scores were remarkably consistent despite varying sequencing depth, contig count, and coverage. ConclusionMetaCompare 2.0 successfully ranked a wide array of environments according to both human health and ecological resistome risks, with both scores being strongly impacted by anthropogenic stress. We packaged the improved pipeline into a publicly-available web service that provides an easy-to-use interface for computing resistome risk scores and visualizing results. The web service is available athttp://metacompare.cs.vt.edu/ 

Abstract Horizontal gene transfer (HGT) occurring within microbiomes is linked to complex environmental and ecological dynamics that are challenging to replicate in controlled settings. Consequently, most extant studies of microbiome HGT are either simplistic experimental settings with tenuous relevance to real microbiomes or correlative studies that assume that HGT potential is a function of the relative abundance of mobile genetic elements (MGEs), the vehicles of HGT. Here we introduce Kairos as a bioinformatic tool deployed in nextflow for detecting HGT events “in situ,” i.e., within a microbiome, through analysis of time-series metagenomic sequencing data. Thein-situframework proposed here leverages available metagenomic data from a longitudinally sampled microbiome to assess whether the chronological occurrence of potential donors, recipients, and putatively transferred regions could plausibly have arisen due to HGT over a range of defined time periods. The centerpiece of the Kairos workflow is a novel competitive read alignment method that enables discernment of even very similar genomic sequences, such as those produced by MGE-associated recombination. A key advantage of Kairos is its reliance on assemblies rather than metagenome assembled genomes (MAGs), which avoids systematic exclusion of accessory genes associated with the binning process. In an example test-case of real world data, use of assemblies directly produced a 264-fold increase in the number of antibiotic resistance genes included in the analysis of HGT compared to analysis of MAGs with MetaCHIP. Further,in silicoevaluation of contig taxonomy was performed to assess the accuracy of classification for both chromosomally- and MGE-derived sequences, indicating a high degree of accuracy even for conjugative plasmids up to the level of class or order. Thus, Kairos enables the analysis of very recent HGT events, making it suitable for studying rapid prokaryotic adaptation in environmental systems without disturbing the ornate ecological dynamics associated with microbiomes. Current versions of the Kairos workflow are available here:https://github.com/clb21565/kairos.